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Abstract 



We consider the scattering of electron by a one-dimensional random poten- 
tial (both passive and active medium) and numerically obtain the probability 
distribution of Wigner delay time (r). We show that in a passive medium 
our probability distribution agrees with the earlier analytical results based 
on random phase approximation. We have extended our study to the strong 
disorder limit, where random phase approximation breaks down. The delay 
time distribution exhibits the long time tail (1/t^) due to resonant states, 
which is independent of the nature of disorder indicating the universality of 
the tail of the delay time distribution. In the presence of coherent absorption 
(active medium) we show that the long time tail is suppressed exponentially 
due to the fact that the particles whose trajectories traverse long distances in 

the medium are absorbed and are unlikely to be reflected. 
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The concept of delay time which was introduced by Wigner way back in 1955 has 
received renewed interest in recent years The delay time in the scattering process 

is generally taken to be related to the duration of collision event or the time spent by the 
particle (the wave packet) in the region of interaction. Very recently universal parametric 
correlations of phase shifts and delay times in quantum chaotic scattering regimes have 
been reported [H. By assuming the validity of random matrices for describing the statistics 
of closed chaotic systems the distribution of inverse delay time is given by the Laguerre 
ensemble from random matrix theory. The delay time statistics is intimately connected 
with the issue of dynamic admittance of microstructures 0. The case study being the 
mesoscopic capacitance [jlO[ fluctuation for the case of a chaotic cavity coupled capacitively 



to the backgate . This has involved calculation of delay time given by the energy derivative 
of the phase associated with the scattering matrix. 

It may be noted that for the case of a single channel the distribution of the delay time for 
a disordered semi-infinite sample has been obtained earlier by using the invariant imbedding 
approach [Q. The stationary distribution Ps{t) for dimensionless delay time r is given by 

= (eW2_i)(i + r2)' 

where A is proportional to the disorder induced localization length and the most probable 
value of r occurs at r^ax = A/2. At long time tail of the above distribution scales as 
The average value of r is logarithmically divergent indicating the possibility of the 
particle traversing the infinite sample before being totally reflected, presumably due to the 



Azbel resonances [|TT|, which makes Landauer's four probe conductance infinite even for 
a finite sample. If the disordered region is semi-infinite, the reflection coefficient will be 
unity, and the complex reflection amplitude will have the form R = e*^'-^-'. If the wave 
packet is incident on the disordered sample it will not be immediately reflected back into 
the lead region, but will be delayed by time proportional to r = hd9/dE. This energy 
dispersive backscattering (or energy dependent random time delay) leads to non-cancellation 
of the instantaneous currents at the surface involving the incident and reflected waves. This 



non-cancellation is expected to lead to a low temperature 1// type noise that should be 
universal (Q-Q. A very recent study claims the delay time distribution in one-channel case 
to be universal and is independent of the nature of disorder (especially the long time tail 
distributions) ^ . We would like to examine this claim through our study. We would like to 
emphasize that to obtain Ps{t) (Eqn. |I|) earlier studies invoke several approximations such 
as random phase approximation (RPA), which is only valid in the small disorder regime and 
moreover, correlation between the phase and the delay time is neglected. 

In view of the recent claim of universality in the delay time distribution, we have studied 
numerically the one-dimensional disordered medium and computed the -Ps(t). Our results 
agree qualitatively with those obtained earlier even in the strong disorder limit where random 
phase approximation breaks down. We also observe same scaling behavior of the tail of the 
stationary delay time distribution for different distribution of the disorder potential. We have 
also studied the delay time distribution in a random active medium, i.e., in the presence of 
coherent absorption. In this case we observe that long time tail distribution is strongly 
suppressed. This is due to the fact that in the absence of absorption the tail distribution 
arises from the possibility of particle traversing the infinite sample (due to resonances) before 
being reflected 0,|TT]]. As these paths traverse a long distance, in the presence of absorption 
the particles will get absorbed in the medium and they cannot be reflected. Consequently 
distribution at the tail gets exponentially suppressed. There are several physical situations 
where one encounters the absorption of elementary particles (electrons, excitons or magnetic 
excitations, etc.) due to impurities or trapping centers in a medium. One recent example 



being the light (photon) propagation in a lossy dielectric medium |[T^-[T4 . 

To describe the motion of a quasi particle on a lattice, we use the following Hamiltonian 
in a tight-binding one-band model: 

H = ^n\n){n\ + V{\n){n + l\ + \n){n-l\), (2) 

where \n) is the non degenerate Wannier orbital at site n, is the site energy at the 
site n and V is the hopping matrix element connecting nearest neighbors separated by 



a unit lattice spacing. We consider three different kinds of disorder where the site ener- 
gies e„ are assumed uncorrelated random variables having distributions which are uniform 
(P(e„) = l/W), Gaussian (-P(e„) oc e"""/^^) and exponential (-P(e„) oc e"""/^). In case of 
modeling absorption we make the site energies complex (e„ = + irj) with the imaginary 
part 7] chosen to be spatially constant (coherent absorption) [Q. The site {n = 1 to N) 
disordered ID sample is embedded in a perfect infinite lattice having all site energies zero. 
The well known transfer-matrix method ||T^ , [l5| . |T6[| is used to calculate the reflection ampli- 
tude r{E) = \r\e~''^^^^ and its phase 9{E) at two values of incident energy E = Eq ± 5E. 
The delay time is then calculated using the definition r = Udd/dE. Throughout our follow- 
ing discussion we consider delay time r in a dimensionless form by multiplying it with V 
and we set = m = 1. 

In view of the fact that the value of incident energy Eq will not change the physics of 
the problem, in the following we choose Eq = Q and dE = 26E = 0.002. In calculating the 
distribution of various quantities we have taken at least 10 000 realizations of the disordered 
sample. The disorder strength (W) and absorption strength [i]) are scaled with respect to V 
i.e., W = W/V and t] = rj/V . To calculate the stationary distribution of delay time we have 
considered a sample of length (L) 5 times the localization length (^), where the localization 



length is calculated by a standard prescription [|T7l- We have confirmed that by changing 
the length L of the sample the distribution remains unchanged. 

In Fig. Ill (a) and (b) we show the plot of numerical data (thin line) for the stationary 
distribution Ps{t) of delay time r for weak disorder {W = 0.5) and strong disorder {W = 2.0) 
respectively. The thick line in the figure is the numerical fit obtained by using the expression 
for Ps{t) given in Eqn. |]. We see that the fit is fairly good even for strong disorder {W = 2.0) 
for which the stationary distribution of the phase of reflected wave, Ps{0), shows (inset of 
Fig. m (b)) two distinct peaks indicating the failure of random phase approximation (RPA) 
in this regime |]18|,|l9[. The values of A thus obtained, when plotted against l/W"^, indicate 
that A scales as l/W"^ as shown in the Fig. ^ Since the energy of incoming particles is fixed, 
the most probable traversal time will be proportional to the typical length traveled by the 



particle in the sample, i.e., to the localization length. As the localization length for passive 
disordered systems |[T^ scales as and A = Tmaxf^, the proportionality of A to l/W"^ 

stands clear. In the inset of Fig. ^ we show the plot of A versus l/W"^ for Gaussian and 
exponential disorder. We see that in the weak disorder limit A scales as l/W"^ for Gaussian 
and exponential disorder as well. 

We would now like to take a closer look at the tail of the stationary distribution Ps{t). 
As mentioned earlier, within RPA the long time delay distribution would scale as The 
numerical least square fit of the expression a/r^ to the long time tail data gives /? ~ 2. This 
is shown in Fig. |^. The appearance of such a tail is attributed to the presence of resonant 
states is due to certain realizations, as in these cases the particle travels a long distance 
before getting reflected . 

Motivated by the recent claim of universality of delay time distribution in one-channel 
case (within RPA), we numerically study the delay distribution Ps{t) for the case of Gaus- 
sian and exponential disorder. The expression in Eqn. |1] was derived within RPA using a 
Gaussian, delta-correlated random potential with zero mean. We have already seen that it 
fits very well with the data for Ps(t) in case of uniform disorder. Similarly, The agreement 
of Eqn. |T] with numerical data for Gaussian and exponential disorder is also excellent in the 
weak disorder limit. 

We now look at the tail of delay distribution and its universality for the three different 
kinds of disorder beyond RPA. Since the origin of tail is due to the appearance of resonant 
realizations which are independent of strength and the type of disorder, we expect that the 
tail distribution would be universal beyond RPA too. In Fig. ^ we plot the tail distribution of 
Ps(r) for uniform, Gaussian and exponential disorder characterized by the strength W = 1.0. 
It can be readily noticed that in all the cases rescaled graphs fall on same curve (within our 
numerical error) clearly indicating the universal nature of tail distribution. For the value 
W = 1.0, we are in a regime beyond RPA as can be seen from the non-uniformity of the 
stationary distribution Ps{6) of the phase of the reflected wave shown in the inset of the 
Fig. ^.Therefore, our numerical simulation results suggest the existence of universality in 
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the long time tail distribution beyond RPA. 

Finally we consider the stationary delay distribution for a coherently absorbing random 
medium. We take the absorption strength rj to be 0.1 and the disorder strength W to be 
1.0 for the active random medium. From Fig. |^ we see that the delay distribution for active 
(coherently absorbing) medium (shown with thick line) falls off much rapidly as compared 
to that of passive medium(shown with thin line). The numerical fit of the tail distribution 
for the active medium, shown in the inset of Fig. ^, indicates that the fall off is exponential 
as Ae~"'^/r^'^. This can be interpreted in terms of suppression of the resonant realizations 
due to the presence of absorption. In case of these resonant realizations, the time spent by 
the particle inside the sample is large as it travels large distance before getting reflected. 
This enhances the probability of the particle getting trapped or absorbed (exponentially) 
in the traversal length thereby suppressing the tail. Analytical study of the distribution of 
delay time in the presence of coherent absorption is under active consideration by us. 

In conclusion, we have numerically studied the universality of the long time tail of the 
stationary delay time distribution for a disordered one dimensional sample. Our study 
reveals that the universality (with respect to the type of disorder) holds for weak as well as 
strong disorder, i.e., irrespective of the validity of the RPA. Also the typical long tail (~ 
1/r^) in the stationary distribution gets exponentially suppressed if we consider a coherently 
absorbing random medium. 
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FIG. 1. The stationary distribution of delay time -Ps(t) for (a) weak disorder (W=0.5) and (b) 
strong disorder (W=2.0). 
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FIG. 2. Plot of A versus l/W"^ for uniform disorder. The numerical fit is valid only for small 
values of disorder strength i.e., W <^ 1.0. The inset shows the plot of A versus for Gaussian 

and exponential disorder. 
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FIG. 3. The plot of tail of Ps(t) for W = 0.5. Numerical fit (thick line) using expression a/r^ 
yields (3 !^ 2.0. 
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FIG. 4. The plot of tail of -Ps(r) for the case of uniform, Gaussian and exponential disorder. 
The disorder strength in all three cases is = 1.0. 
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FIG. 5. The plot of -Ps(t) for passive {rj = 0, shown with thin hne) and active (r/ = 0.1, shown 
with thick hne) disordered medium. The inset shows the tail of Ps{t), numerical data using thin 
line and numerical fit using thick line, for the absorbing medium. 
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